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ABSTRACT 

The distribution of QSO radio luminosities has long been debated in the literature. Some 
argue that it is a bimodal distribution, implying that there are two separate QSO populations 
(normally referred to as 'radio-loud' and 'radio-quiet'), while others claim it forms a more con- 
tinuous distribution characteristic of a single population. We use deep observations at 20 GHz 
to investigate whether the distribution is bimodal at high radio frequencies. Carrying out this 
study at high radio frequencies has an advantage over previous studies as the radio emission 
comes predominantly from the core of the AGN, hence probes the most recent activity. Studies 
carried out at lower frequencies are dominated by the large scale lobes where the emission is built 
up over longer timescales (10 7 — 10 8 yrs), thereby confusing the sample. Our sample comprises 
874 X-ray selected QSOs that were observed as part of the 6dF Galaxy Survey. Of these, 40% 
were detected down to a 3er detection limit of 0.2-0.5 mJy. 

No evidence of bimodality is seen in cither the 20 GHz luminosity distribution or in the 
distribution of the R20 parameter: the ratio of the radio to optical luminosities traditionally used 
to classify objects as being either radio-loud or radio-quiet. Previous results have claimed that 
at low radio luminosities, star formation processes can dominate the radio emission observed in 
QSOs. We attempt to investigate these claims by stacking the undetected sources at 20 GHz and 
discuss the limitations in carrying out this analysis. However, if the radio emission was solely due 
to star formation processes, we calculate that this corresponds to star formation rates ranging 
from ~ 10M Q yr" 1 to ~ 2300M© yr" 1 . 

Subject headings: quasars: general — galaxies: active — galaxies: star formation — radio continuum: 
galaxies 

1. Introduction 

1 Present address: ASTRON, the Netherlands Institute 

for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The QSOs are often classified into two broad cate- 

Netherlands. gories based on their radio properties: radio-loud 

and radio-quiet, where the 'radio-loudness' of a 
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QSO is usually defin ed by the ratio of its r adio 
to optical luminosity (jKellermann et al . 1989 and 
references therein). However, the underlying dis- 
tribution of QSO radio luminosities has long been 
debated in the literature. There are two opposing 
views: the first is that the distribution is bimodal 
with approximately 5-10% of QSOs bei ng radio- 
loud (jlvezic et al.ll2002l : I Jiang et al.ll2007t ). and the 
second is that there is a broad, continuous distri- 
bution with no clear dividi ng line between radio- 
loud and radio-qu iet QSOs ( Cirasuolo et al. 20031: 
Rafter et al.ll2009h . 

Clarifying this issue can provide insight into the 
physics associated with forming radio jets. A bi- 
modal distribution suggests that there are two in- 
trinsically different classes of QSOs, only one of 
which is able to become a strong radio source, 
while a continuous distribution suggests that all 
QSOs have low-luminosity radio sources which be- 
come stronger during episodes of unusually high 
activity. 

Although there have been numerous stud- 
ies addressing this issue, the wide range of se- 
lection criteria used to define the samples ob- 
served, along with the different flux density 
limits reached, contribute to a large range of 
contradicting results. Studies that report a 
bimodal distribution (i.e. two differe nt QSO 



Stocke et al 



populations) incl ude Strittmatter et al] ( 1980h 



Kellermann et all (fl989h: iMiller et all 



1990); 



(119921): iHooper et all (jl995l ) and 
Ivezic et alj ~ ll2002t l2004h . On the other hand 



studie s bv ICondon et al. (Il980l) : IWhite et al 



2000): Cirasuolo et alJ 



(2003) 



and iRafter et al 



2009) argue that the range of radio luminosities 



produces a more continuous distribution which in- 
cludes a population of 'radio-intermediate' QSOs. 

With the completion of large area surveys 
such as the Sloan Digital Sky Survey (SDSS; 



York et all l2000h . Faint Images o f the Radio 



Sky a t Twenty centimetres (FIRST; I Becker et al.l 
19951) and the NRA O VLA Sky Survey (NVSS; 
Condon et al. 19981) came a shift from follow-up 
observations of smaller samp les to crossmatch- 
ing of large catalog s (e.g 



Ivezic et al 
20071) . 



2002; 



Cirasuolo et al.1 120031: Uiang et al.l 120071) . These 
studies have led to suggestions that the frac- 
tion of radio-loud sources observed is correlated 
with a range of observed properti es including the 
optical magnitude and redshift (jPadovanil 11993 ; 



Jiang et al] 120071 ). optic al colors 



20071 ). black hole mass (IBest et al 
cretion rate { Sikora et alF 20071 



White et al 
20051) 



ac- 



Shankar ct al 



20101: iBroderick fe Fenderll201 ll) . host galaxy mor 



phology (IStawarz et al.l 120081) and ste llar mass 
(IBest et al.ll2005t iMauch fc Sadlerll2007l) . 



Another reason why different studies have pro- 
duced contradictory results is that many of the 
previous radio studies have used QSO samples 
which span a wide range in redshift and a rela- 
tively narrow range in optical magnitude. This 
introduces a strong and spurious correlation be- 
tween redshift and optical luminosity, which in 
combination with the rapid cosmic evolution of the 
QSO population introduces selection effects which 
greatly complicate any interpretation of the data 
(| Jiang et alll2007l ). 

By carrying out this study at high radio fre- 
quencies, we probe the radio emission much closer 
to the core of the AGN. As such, we can gain in- 
sight into the most recent activity, comparable to 
timescales seen in the optical regime. At lower 
radio frequencies the large scale radio lobes domi- 
nate the emission which are built up on timescales 
of 10 7 — 10 8 yrs and do not necessarily relate to 
the current AGN activity. 

Upgrades of existing radio facilities to signifi- 
cantly increased bandwidth such as the Compact 
Array Broadband Backend (CABB; IWilson et al 



201 lh on the Australia Telescope Compact Ar- 



ray (ATCA) and the Expand ed Very Large Ar- 
ray (EVLA: IPerlev et al1l201ll ) has enabled stud- 
ies of large numbers of objects to much deeper 
flux density limits, probing further into the radio- 
quiet regime. This has recently been shown in 
Kimball et al] ( 2011 ) who observed a sample of 



179 optically-selected QSOs in the redshift range 
0.2 < z < 0.3 with the EVLA down to an rms 
of 6-8^iJy. The authors detect nearly all sources 
in the sample (97%) and find that the resulting 
radio luminosity function is consistent with 2 ra- 
dio source populations: those dominated by AGN 
emission and those dominated by star formation 
in the host galaxy. 

In this paper, we observe a sample of X-ray se- 
lected QSOs with z < 1 that were obse rved as part 
of th e 6dF Galaxy Survey (6dFGS; Ijones et aT 



2004) at 20 GHz to study the bimodality of the 



radio luminosity distribution of X-ray selected 
QSOs at high radio frequencies. Section [2] de- 
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tails the sample selection and Section [3] describes 
the follow-up radio observations. Results of this 
analysis are presented in Section |4l followed by a 
discussion on stacking the undetected sources in 
Section [5j We conclude in Section [B] 

2. Sample Selection 

Previous studies addressing the dichotomy of 
radio-loud and radio-quiet QSOs have generally 
been based on optically selected samples which are 
then followed up at radio wavelengths. We take a 
slightly different approach by starting with an X- 
ray selected AGN sample. Target s were selected 



from the RASS-6dFGS catalog (jMahonv et al 
2010t ); a catalog of 3405 AGN selected from 



the ROSAT A LL Sky Survey (RASS) Bright 
Source catalog ( Voges et al. 19991) that were ob- 
served as part of the 6dF Galaxy Survey (6dFGS; 



Jones et al.l 120041 ) . The 6dF Galaxy Survey is a 



spectroscopic redshift survey covering the entire 
southern sky with \b\ > 10°. Bright X-ray sources 
were selected such that AGN activity was the pri- 
mary source of X-ray emission, rather than X-ray 
emission from clusters. Of these, 2224 were ob- 
served as part of the 6dFGS3 of which 77% have 
reliable spectroscopic redshifts. For sources where 
a reliable spectrum was not able to be obtained, 
this was primarily due to low signal-to-noise ra- 
tios such that absorption features could not be 
identified. These sources would not be included in 
our QSO sample due to the lack of broad emission 
lines (as described below) so the 23% of sources 
without reliable redshifts should not add any sig- 
nificant additional bias. 

To ensure the sample comprised solely QSOs, 
only targets with M&j < —22.3 and broad emis- 
sion lines in the optical spectrum were selected^. 
This absolute magnitude cutoff ensures that the 
optical luminosity is dominated by nuclear emis- 
sion with negligible contribution from the galaxy 
continuum. We also placed a redshift cut of z = 1 
which leaves us with a manageable sample, yet still 



1 This X-ray selected sample was an additional target sample 
in the 6dF survey meaning that sources were observed with 
any spare fibres that were available 

2 This absolute magnitude limit corresponds to the stan- 
dard QSO crite ria of Mjj < —23 introduced by 

I Schmi dt & GrceiJ ||1983T ) shifted to current cosmological 
parameters with Hq = 71kms~ 1 Mpc - 1 and Q m = 0.27. 



allows us to investigate any trends with redshift. 
The final sample contains 874 objects. 

3. High-frequency Radio Observations 

In contrast to previous studies which gener- 
ally use radio emission at either 5 or f.4GHz to 
determine the radio luminosity distribution, we 
have observed our sample at a higher frequency of 
20 GHz. Carrying out this study at 20 GHz means 
that the core of the AGN is the dominant source 
of emission and provides insight into the most re- 
cent AGN activity. At lower radio frequencies, the 
emissions observed is dominated by the radio lobes 
which could be relics of past activity, thereby con- 
fusing the sample. 

3.1. The AT20G Survey 

Before reobserving the 874 targets which com- 
prise the full sample, we crossmatched the sam- 
ple with the Australia Telesc ope 20 GHz survey 
(AT20G; iMurphv et all l201dh to search for pre- 
vious detections. The AT20G survey is a blind 
survey of the southern sky down to a flux density 
limit of 40 mJy. The final catalog consists of 5890 
sources and is domina ted by flat spectrum sources 
(jMassardi et al.ll201lf) . We found 56 sources in our 
X-ray selected sample that were detected in the 
AT20G survey and therefore were not followed up 
any further. 

3.2. New Observations 

The remaining 818 X-ray selected QSOs were 
observed with the Australia Telescope Compact 
Array (ATCA) from 2008-2010. The observations 
were carried out using the compact, hybrid config- 
uration (H168) over four separate observing runs. 
These are summarised in Table [T] 

The first observing run carried out in October 
2008 used the old correlator that had a band- 
width of 128 MHz in each of the two IFs (here- 
after denoted as 2x128). Central frequencies of 
18752 and 21056 MHz were chosen in order to 
match the observing setup of the AT20G survey. 
In early 2009 the new Compact Array Broadband 
Backend (CABB) was installed on t he ATCA, in- 



creas ing the bandwidth to 2 x 2 GHz (jWilson et al 



20111 ). The remaining three observing runs used 
this new setup with central frequencies of 19000 
and 21000 MHz. 
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The observations were carried out in a two-step 
process; all sources were observed for 2 x 40 sec cuts 
and targets not detected in this time were then 
reobserved for 2x5 min. A small area of sky was 
affected by poor weather in all scheduled observing 
runs and therefore 18 sources were unable to be 
observed. These sources have been removed from 
the sample. 

3.3. The deep sample 

Due to poor weather during the observations we 
were unable to observe all sources for the 2x5m 
integration time intended. A 'deep sample' was 
defined by a region of sky (shown in Figure [1]) in 
which sources that were not detected in 2x40s 
were reobserved for an additional 2 x 5 m. The 
deep sample covers most of the sky south of —30° 
declination and a few other regions which were ob- 
served during good weather. Note that the deep 
sample is defined only by the area of sky covered 
and includes sources that were detected in the ini- 
tial 2x40s observations. This ensures there is no 
bias by only including the fainter radio sources. 
There are 397 objects in the deep sample. 

3.4. Data reduction 



All data were reduced using MIRIAD ( Sau lt et al 



1995). Images were deconvolved only if a 3a peak 
was detected at the optical position of the source, 
or a 5cr peak elsewhere in the field. Since we 
are only interested in emission from the core of 
the AGN, of which the position is known from 
the 6dF observations, fluxes and associated errors 
were obtained by fitting point sources at the op- 
tical position in the deconvolved images using the 
MIRIAD task IMFIT. As the optical position was 
known to better than 1" accuracy, a detection at 
20 GHz was accepted if the flux density exceeded 
the 3cr threshold. 

In cases where other radio sources were detected 
in the field, these sources were modelled in the 
(u,v) plane and removed from the visibilities. The 
field was then imaged again with these sources re- 
moved. This removed any sidelobes which may 
affect the flux density measurement of the central 
source. Due to the small field of view and low 
source density at 20 GHz, the majority of images 
had at most one other radio source in the field that 
needed to be removed. In a very small number of 



cases this process meant that radio lobes possi- 
bly associated with the X-ray /optical source were 
also removed. However, as we are only interested 
in the core emission to probe the most recent ac- 
tivity as opposed to emission from previous epochs 
of activity, this does not affect the results of this 
study. 

4. Results 

4.1. Detection rates 

Of the 800 sources observed, 321 (40%) were 
detected above the 3a threshold. Figure [2] shows 
the flux density distribution of both the full sam- 
ple and the deep sample. Detections are shown by 
the solid line and 3a upper-limits shown by the 
dashed line. The distribution of the upper-limits 
in the full sample clearly shows the effect of the 
different integration times for the sample. The 
non-detection peak at approximately 0.2 mJy cor- 
responds to those sources that were observed for 
the extra 2x5 minute integrations and the peak at 
~0.5mJy corresponds to those sources that were 
only observed for 2x40 second integrations. For 
the 397 objects in the deep sample, 187 (47%) were 
detected. 

Redshifts were known for all sources in the sam- 
ple and radio luminosities were calculated. We 
assumed a flat spectral index of a — when 
applying K-corrections since we are focussing on 
the 20 GHz radio emission that comes predomi- 
nately from the core of the AGN. The median 
spectral indices of sources in the AT20G survey 
are a\ = -0.16, a% = -0.16 and af = -0.28 
indicating that high frequency radio sources have 
flat spectral indices o ver a large frequency range 
(jMassardi et al]|2011l) . Assuming a spectral index 



of 



,20 



-0.5 instead, only shifts the luminosity 



by 0.02 (in logspace) for a source at z = 0.1 and 
0.15 for a source at z = 1. The 20 GHz luminosi- 
ties as a function of redshift are shown in Figure 

El 

As with any flux density limited sample, this 
sample is subject to a luminosity-redshift degen- 
eracy as can be seen in Figure [3] To better illus- 
trate the selection effects of the sample, Figures 
|4] and [5] show the X-ray and optical luminosities 
respectively as a function of redshift. 

To quantify the 'radio-loudness' of an object we 
have also calculated the R20 parameter; the ratio 
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Table 1: List of observations for this program. All the observations from 2008-2010 used the Hybrid 168m 
array configuration to obtain better (u,v) coverage. The October 2008 run was using the old correlator 
(2x128 MHz bandwidth) and all other runs used CABB (2x2GHz bandwidth). 



Date of 


Central 


Bandwidth 


Time on 


approx. 3<T 


Observations 


Frequencies (MHz) 


(MHz) 


Source 


Detection limit 


October 2008 


18752, 21056 


2x128 


2 x40s 


3 mjy 


April 2009 


19000, 21000 


2x2000 


2 x40s 


0.5 mjy 


October 2009 


19000, 21000 


2x2000 


2 x40s 


0.5 mjy 




19000, 21000 


2x2000 


2 x5m 


0.2 mjy 


March 2010 


19000, 21000 


2x2000 


2 x5m 


0.2 mjy 




Fig. 1. — Polar projection of the source distribution. The black squares indicated sources that were detected 
in the AT20G sample and therefore were not followed up as part of these observations. Green triangles were 
sources that were not observed due to poor weather and therefore removed from the sample. The blue circles 
are objects in the 'deep sample' which contains sources that, if not detected in 2x40 s, were reobserved for 
an additional 2x5m. The area of sky defining the deep sample is shown by the solid line. Red diamonds 
indicate those sources that were only observed for 2x40s and are included in the 'full sample'. 
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140 




full sample. The blue circles are sources which 
were detected at 20 GHz, while the red triangles 
Fig. 2. — Flux density distribution for the full indicate 3a upper limits, 

sample (top panel) and the deep sample (bottom 
panel), binned in log flux. Sources that were not 
detected arc shown as upper limits by the dashed 
histogram and detections shown by the solid his- 
togram. 
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Fig. 3. — 20 GHz radio luminosity for the full sam- 
ple as a function of redshift. The blue circles are 
sources which were detected at 20 GHz, while the 
red triangles indicate 3<r upper limits. 




Fig. 5. — Optical luminosity against redshift for 
the full sample. Again the blue circles are sources 
which were detected at 20 GHz, while the red tri- 
angles indicate 3<r upper limits. 
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of the radio to optical luminosities: 



-R20 = log(y^) 



(i) 



The bj-band magnitudes were used to calculate 
the optical luminosities, and shifted into the rest- 
frame using a spectral index of a — —0.5, typical 
for QSOs. Throughout this paper we use the fol- 
lowing cosmological parameters: H = 71 km s -1 
Mpc" 1 , Sl m = 0.27 and Q\ = 0.73 (jLarson et al 
20111) . 



4.2. Radio luminosity distributions 

The 20 GHz radio luminosity distributions for 
both the full sample (top panel) and deep sample 
(bottom panel) are shown in Figure [6] Sources 
not detected are shown as 3er upper limits by the 
dashed line. The ambiguity surrounding whether 
the luminosity distribution is bimodal or not also 
means that there are varying definitions for classi- 
fying an object as 'radio-loud' in the literature. An 
often used dividing line separating radio-loud and 
radi o-quiet sources is defined as P2n = 10 24 W/H z 
(e.g. iMiller et alJll990HGoldschmidt et alll9 99il 
Using this luminosity as a approximation to sepa- 
rate the radio-loud and radio-quiet sources means 
that although not all sources were detected, our 
observations were sensitive enough such that all 
of the sources that weren't detected fall into the 
radio-quiet regime. As such, any minimum in the 
distribution corresponding to the division of radio- 
loud and radio-quiet sources would be seen. Figure 
[6] shows no indication of a bimodality. 

The distribution of the R20 parameter for 
both the full and deep samples is shown in Fig- 
ure [7] with upper limits again shown by the 
dashed line. Here the division between radio- 
loud and r adio-quiet sources is typ i cally given as 
R20 = 1 (Kellermann et al. 19891 : Hooper et al 



1995t I Jiang et al.l 120071 ) . again meaning that we 
have detected all the radio- loud sources. Figure [7J 
also shows no clear evidence for a bimodal distri- 
bution. 



3 The division at 10 24 is usually defined at frequencies of 1.4 
or 5 GHz, but using the assumption that QSOs have a flat 
spectrum, it will be the same at 20 GHz 




log 20 GHz radio luminosity (W/Hz) 



Fig. 6. — Radio luminosity distribution for all 
sources (top panel) and the deep sample (bottom 
panel) . 





i_ _ 2 upper limits 




1 1 full sample 








R 2 „ parameter (log(Lj„/L„„)) 

Fig. 7. — R20 distribution for all sources (top 
panel) and the deep sample (bottom panel). The 
R20 parameter is defined as the log of the ratio of 
the 20 GHz rest frame luminosity and the optical 
(rest frame bj-band) luminosity. 
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4-2.1. Selection effects 



One major limitation in compiling a flux den- 
sity limited sample that covers a large range in 
redshift is that we are sensitive to different lumi- 
nosities as a function of redshift. At high z we 
only detect the most luminous sources while at 
low z we observe the more common, less luminous 
sources, but miss the bright sources due to the 
small volume sampled. 

This effect is clearly shown in Figure [8] which 
plots the 20 GHz radio luminosity distribution 
split into separate redshift bins. In every redshift 
range there is a broad range of radio luminosi- 
ties observed, but no evidence for any bimodality. 
Figure [9] shows the R20 distribution split into sep- 
arate redshift bins where we see the same redshift- 
luminosity degeneracy. There is a clear shift from 
a predominately radio-quiet population at low red- 
shifts, to a radio-loud population at high redshifts 
due purely to the flux density limit of the sample. 

Although we are using an X-ray selected sam- 
ple, applying the condition that all objects have 
a 6dF spectrum adds an additional optical flux 
limit. This optical flux limit is made even stricter 
by adding the condition that the absolute mag- 
nitude must be brighter than M&j = —22.3 lead- 
ing to a narrow range of optical luminosities ob- 
served. While the sample presented here covers 
a relatively large redshift range from < z < 1, 
many previous studies place no redshift limit at 
all. It is this redshift-luminosity degeneracy, plus 
the numerous different selection criteria, that have 
most likely lead to contradictory results in previ- 
ous studies on this subject. 

4.2.2. The validity of the R parameter 

It has previously been shown that the R param- 
eter is a meaningful measure of the radio loudness 
only if the optical and radio lum inosities are corre- 
lated (|Goldschmidt et al.|[l999l) . Figure ITUl shows 
that this is not the case. Our sample covers a wide 
range in radio luminosity (more than 8 orders of 
magnitude), but a much narrower range in optical 
luminosity (~2 orders of magnitude). 

However, since the R parameter is the ratio of 
optical to radio luminosities, it does provide useful 
information on the underlying physics operating in 
these systems. Taking the optical luminosity as a 
proxy for the accretion power and the radio lumi- 
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Fig. 8. — 20 GHz radio luminosity distribution for 
the full sample separated into redshift bins. There 
is no evidence for a bimodal distribution in any 
redshift range up to z = 1. 




parameter (logfLjj/L )) 



Fig. 9. — R20 distribution for the full sample sep- 
arated into redshift bins. 
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nosity as a proxy for the jet power, from the ratio 
we can gain insight into the accretion efficiency 
of these objects. The fact that we observe a nar- 
row range in optical luminosities, but a very broad 
range of radio luminosities indicates that the ac- 
cretion power does not map directly to the radio 
jet power observed. This suggests that other pa- 
rameters (such as the magnetic field or the spin 
of the central black hole) can play a large role in 
determining the radio output observed in QSOs. 

5. Discussion 

The continuous distribution of high-frequency 
radio luminosities suggests that radio-loud and 
radio-quiet QSOs are not two intrinsically differ- 
ent populations. However, it has recently been 
suggested that while the distribution of radio lu- 
minosity of QSOs does not appear bimodal, it can 
be fitte d by two populations : AGN and star for- 
mation ( Kimball et al.l 12011 ). This was based on 
a study of an optically selected sample of QSOs 
that were observed at 6 GHz using the EVLA. 

Unfortunately the large number of non de- 
tections in our sample makes it difficult to di- 
rectly compare our 20 GHz luminosi ty distribution 



with t he distribution presented in Kimball et al 
( 201ll ). To investigate this further requires more 



information on those sources below our detection 
limit. 

5.1. The limitations of stacking 

One of the most obvious ways to gain this extra 
information on the undetected sources is to stack 
all of the non detections at the position of the 
QSO to obtain an average flux density. To ensure 
that the noise is approximately the same in each 
of the individual observations, only objects in the 
deep sample were stacked. This resulted in a 16er 
detection of 0.074±0.005mJy. 

However, rather than being an indication of the 
underlying flux density of the sample, this stacked 
emission is strongly biased towards those sources 
just below our detection limit. As such, this stack- 
ing detection is more indicative of our survey flux 
density limit rather than being an intrinsic prop- 
erty of the undetected QSOs. In addition, due 
to the large range in redshifts combined into the 
one image, further complications arise if we try to 
convert this measured average flux density into a 



luminosity. The large range in redshift also means 
that each image is weighted differently and simply 
taking the average or median luminosity distance 
will not have the same weighting. 

The very nature of stacking means that we lose 
all information on the underlying distribution of 
the non detected sources and therefore it is not 
trivial to compare the stacked dete ction with the 



lumin osity distribution presented in lKimball et al 
(1201 lh . 



An alternative method is to plot the flux den- 
sity probability distribution function (hereafter 
termed flux density PDF) as shown in Figure [TT] 
This is derived by measuring the flux densities at 
the optical positions of the targets, shown by the 
solid line. Sources with more than a 3tr detec- 
tion at 20 GHz are shaded in grey. We have again 
only included sources in the deep sample as these 
images have approximately uniform noise proper- 
ties. Although the non-detections do not have a 
significant flux density measurement in the indi- 
vidual images, this flux density PDF shows that 
the average flux density of the undetected QSOs 
is non-zerc@. 

To confirm that this is not a result of the noise 
in the images, we have also calculated the distri- 
bution taken from random points in the radio im- 
ages, shown by the dashed line. The median flux 
density of the PDF for our X-ray selected QSOs is 
0.165 mJy, compared to the median of -0.001 mJy 
for the random distribution. The median flux den- 
sity of the undetected sources is 0.074 mJy, iden- 
tical to the flux density obtained from the stacked 
image. 

While it is again non-trivial to convert this dis- 
tribution into luminosity space, this does provide 
information on the underlying distribution as op- 
posed to simply obtaining an average flux density. 
With some sophisticated modelling it would be 
possible t o take the QSO lumin osity function pre- 
sented in Kimball et al. I (|201lh and predict what 
would be observed at 20 GHz by assuming either 
a flat or steep spectrum population and compare 
to the results presented here. However, this is be- 
yond the scope of this paper and will be addressed 



4 When carrying out the radio observations, the optical po- 
sition of the QSO was slightly offset from the phase center, 
so this non-zero peak in the distribution is not due to a DC 
offset at the phase center. 
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in future works. 



v upper limits 
• detections 




log optical luminosity (W/Hz) 



Fig. 10. — 20 GHz radio luminosity against bj- 
band optical luminosity. The dotted lines mark 
lines of constant R20. 
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Fig. 11. — Statistical probability distribution for 
our X-ray selected sample (solid line) which shows 
that while many sources are not detected in the in- 
dividual images the average 20 GHz flux density of 
the sample is non-zero. The dashed line shows the 
distribution calculated from random points which 
confirms that the QSO distribution is real, and 
the dotted line denotes the median flux density of 
0.165 mjy. Only sources with a peak flux density 
less than 2 mjy are shown here, but the tail of the 
distribution extends to brighter flux densities. 



5.2. Star formation rates 

Although we cannot distinguish whether the ra- 
dio emission is dominated by AGN or star for- 
mation processes in the data presented here, we 
can calculate rough estimates of the star formation 
rate (SFR) expected using the following relation: 



SFR(L 1A ) 



Li. 



(2) 



8.85 x 10 20 

where the SFR is in M© yr and L1.4 is the 
1.4 GHz radio luminosity in W/Hz. This relation 
assumes a constant burst of star-formation over 
100 Myr dSullivan et al.ll200lh . 

Assuming a steep spectral index of a = —0.7 
and using a 20 GHz flux density of 0.074 mjy (the 
average flux density of the undetected sources) , we 
can extrapolate a 1.4 GHz flux density of 0.48 mjy. 
At a redshift of z = 0.1 this gives a rough SFR 
of ~ 10M Q yr -1 and at z = 1 this would be 
~ 2300Mq yr -1 . Note that these estimates as- 
sume that all of the radio emission is due to star 
formation processes and so only represent upper 
limits of the star formation rates. 

5.3. QSO evolution models 

In current galaxy evolution models, QSO ac- 
tivity is attributed to the accretion of cold-gas 
(|Croton et alJl2006t [Bower et alJl2006h . As such, 
it is reasonable to assume that this same reser- 
voir of gas also contr ibutes to star forma tion. 
Using Herschel data, Bonfield et al.l (|201lh re- 
cently found that the star formation (calculated 
using the Infra-Red (IR) luminosity) is correlated 
with the QSO accretion luminosity, indicating that 
both processes are supplied by the same cold gas. 
This a g rees with CO observat ions o f lScoville et al 



(|2003D : iBertram et al.l (|2007D who find that the 
majority of low luminosity QSOs resi de in gas-rich 



galaxi es. In addition, recent results of lPierce et al 



(|201lh find SFRs ranging from 3-29 Mq yr" 1 for 
a sample of X-ray selected AGN (not necessarily 
all QSOs) providing further evidence that star- 
formation processes play an important role at 
lower radio luminosities. 

According to the current picture of galaxy 
formation through hierarchical merging, gas-rich 
galaxies merge to form Ultra-Luminous Infra-Red 
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Galaxies (ULIRGs). These mergers provide the 
large reservoirs of cold gas needed to fuel the 
AGN (and star- formation) . At this early stage 
the galaxy is shrouded in dust, hence is very lu- 
minous in the IR. As the material begins to settle 
towards the center of the galaxy, AG N activity can 



be switched on, resulting in a QS O (jSanders et al 



19881 : ICanalizo fc Stockton] l200lh . As such, it is 
not unreasonable to expect star formation pro- 
cesses to start to play a dominate role in the radio 
emission observed at lower luminosities. In fact, 
the star formation rate could also provide insight 
into the evolutionary stage of the QSO. During 
the ULIRG phase high star formation rates are 
expected, while lower star formation rates could 
indicate that the gas reservoir is close to being 
depleted. 

5.4. Future Observations 

From the radio data presented here alone we 
cannot distinguish whether the radio emission 
in the low-luminosity QSOs comes from AGN or 
star-forming processes. However, there is mount- 
ing evidence that star formation in the QSO 
host becomes the dominant source of radio emis- 



sion at these luminosities (Kimball ct al.l 12011 
Padovani et~aLll201lb iBonfield et al.ll201lh . 



To properly distinguish whether the emission is 
AGN or star formation dominated, better spatial 
information is needed, ideally using Very Long 
Baseline Interferometry (VLBI). However, it is 
most likely that both processes contribute to the 
observed flux density at some level. As such, 
molecular line observations with the Atacama 
Large Millimeter Array (ALMA) can study the 
distribution of cold gas, potentially disentangling 
the gas being accreted onto the central black hole 
from gas responsible for star formation in the host 
galaxy. 

6. Conclusions 

We have observed a sample of 874 X-ray se- 
lected QSOs that were observed as part of the 
6dF Galaxy Survey at 20 GHz to determine if 
the high-frequency radio luminosity distribution 
is bimodal. High-frequency observations provide a 
unique dataset where the emission is dominated by 
the core of the AGN, hence only showing the most 
recent activity. Deep observations at 20 GHz show 



no evidence for a bimodal distribution in either the 
radio luminosity distribution or in the distribution 
of the R20 parameter (the ratio of optical to radio 
luminosities). The broad range of radio luminosi- 
ties observed for a very narrow range in optical 
luminosity indicates that the accretion power (ap- 
proximated by the optical luminosity) can result 
in a very wide range of jet luminosities. 

Investigating the redshift effects of these distri- 
butions clearly show a strong redshift-luminosity 
degeneracy, which if not taken into account, could 
lead to an apparent bimodal distribution. This ef- 
fect could explain the contradictory results seen in 
previous studies. 

Recent results have suggested that at low lu- 
minosities star formation can dominate the ra- 
dio emission observed in QSOs. To investigate 
this further we attempted to stack the undetected 
sources at 20 GHz. This results in a 16a detection 
of 0.074±0.005mJy, but is heavily biased towards 
those sources just below the detection threshold. 
Instead we derive the flux density probability dis- 
tribution by calculating the peak flux density at 
the position of all sources in the sample, which 
has a median flux density of 0.165 mJy. This pro- 
vides more information on the underlying distri- 
bution than is achieved by stacking. By mod- 
elling the QSO lumi nosity functions presented in 
Kimball et al.l ( 201lh it would be possible to pre- 
dict the 20 GHz luminosity function assuming ei- 
ther a flat or steep spectrum population and com- 
pare with the results presented in this paper. 

While this modelling is beyond the scope of this 
paper, we can obtain rough estimates of the star 
formation rate assuming all of the radio emission 
in the low luminostiy QSOs is due to star forma- 
tion. Using the average 20 GHz flux density calcu- 
lated for the undetected sources, this corresponds 
to star formation rates of ~ 10M Q yr _1 at z = 0.1 
to ~ 23OOM yr" 1 at z = 1. 
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